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Abstract 

We consider the time-dependence of the Higgs vacuum expectation value [vev) given 
>^ _ by the dynamics of the Standard Model and study the non-adiabatic production of 

ly-^ ' both bosons and fermions, which is intrinsically non-perturbative. In the Hartree ap- 

, proximation, we analyse the general expressions that describe the dissipative dynamics 

I due to the back-reaction of the produced particles. In particular, we solve numerically 

■ some relevant cases for the Standard Model phenomenology in the regime of relatively 

. small oscillations of the Higgs vev. 

o 



^ , 1 Introduction 

^ • In the Standard Model, several fundamental constants such as the Fermi coupling G and 
the mass of gauge bosons and fermions depend on the vacuum expectation value (vev) of 
^ ! the Higgs field, because of the well-known mechanism of spontaneous symmetry breaking. 
^ I The equation of motion of the Higgs field on the other hand allows for (periodic) time- 
dependent solutions for the Higgs vev, which can then be viewed as a particular case of 
varying fundamental constants (for a review, see Ref. p[]). This behaviour differs from the 
more common cases of adiabatic variations of the fundamental constants and, due to its 
periodic nature, can lead to efficient (non-adiabatic) particle production. 

The issue of the constancy of the constants of physics was (probably) first addressed by 
P. A. M. Dirac [21 [BJ with his "Large Numbers hypothesis" : very large (or small) dimen- 
sionless universal constants cannot be pure mathematical numbers and must not occur in 
the basic laws of physics. He then proposed that they be considered as (typically slowly) 



^email: casadio@bo.infn.it 
^email: iafelice@bo.infn.it 
email: vacca@bo.infn.it 



1 



variable parameters characterising the state of the Universe and pointed out the possibihty 
of measuring astrophysical quantities to settle this question [H El E] • More recently, theo- 
ries such as string theory and models with extra spatial dimensions have also predicted the 
time-dependence of the phenomenological constants of the low energy regime describing our 
Universe [3 E] . 

The cases we shall consider instead correspond to oscillations of the Higgs vev with 
periods (set by the Higgs mass scale) of the order of 10^^^ s. Such a behaviour can be 
obtained from the usual dynamics of the spatially homogeneous Higgs field [5]. We therefore 
begin by considering the classical equation of motion for the time-dependent Higgs vev (i.e. 
a classical condensate) in a homogeneous patch of space-time and identifying the relevant 
regimes. Quantum fluctuations of fermion and boson fields of the Standard Model are then 
analysed on this Higgs background and explicit expressions for the particle production [TU], 
an intrinsically non-perturbative effect, are presented. In particular, we investigate which 
bosons and fermions are produced more abundantly depending on the Higgs mass which we 
consider in a physically sound range of values [iTj. Further, the back-reaction of particle 
production [12] is analysed in the Hartree approximation which is well suited to describe 
such a dissipative effect. Similar methods have been previously employed to study particle 
production in strong fields [13] and pre- as well as re- heating in Cosmology [TU [THl US]- 

According to our findings, particle production induced by the oscillating Higgs can be 
very efficient. From the phenomenological point of view, if the Higgs were oscillating now, 
the amplitudes of such oscillations should therefore be extremely small. On the other hand, 
this mechanism could explain how Higgs kinetic and potential energy have dissipated in 
the past [Ej. Note that we typically consider regimes such that only small oscillations 
around an absolute minimum are present and the symmetry breaking phase transition is 
not significantly affected. However, since we solve the complete system of coupled non-linear 
equations that describe particle production and their back-reaction, one does not expect that 
the produced particles are thermal [T7]. They also turn out to have mostly small momenta 
and are therefore non-relativistic. 

In Section [21 we review the solutions of the classical equation of motion for the Higgs vev 
to properly identify the periodicity properties. Fermion modes on such a Higgs background 
are then studied in Section [3], with a particular emphasis on the introduction of physical 
quantities and the regions which lead to quantum particle production. The same kind of 
analysis is then given for the bosons in Section jH Fermion and boson production within 
the Standard Model are then studied in details in Sections and [HI respectively, taking into 
account the dependence on the Higgs mass. At this stage we neglect the back-reaction. This 
latter effect, which has to be included to describe consistently the dynamics is then considered 
in Section [71 In this section we analyze phenomena like the dissipation of the Higgs energy 
and particle production, consistent with total energy conservation. Final comments are given 
in Section [HI 

We shall use natural units with h = c = 1. 
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2 Higgs vev time dependence 



The Higgs sector of the Standard Model Lagrangian we are interested in is given by [T8] : 

Ch = d^H^d^H - ^2 H^H - i A {H^Hf, (2.1) 
where H is the complex iso-doublet 

H=^exp (i^] [v + h{x)] (^] , (2.2) 



72 "V V J'^ ' '^^^^^Vl 

t" (with a = 1,2,3) are the SU{2)l generators, C and h (with (0|dO) = = (0|/i|0)) are 
scalar fields which parameterise the fluctuations of H around the vacuum state. We have 
neglected the couplings with gauge and fermion flelds. 

We flnd it useful to introduce two new quantities, Mh and /?, deflned as 

MH = 2—^/X, ^^ = 13 -I Ml, (2.3) 
9 ^ 

where g is the SU{2)l coupling constant and Mh is the "bare" Higgs boson mass (corre- 
sponding to constant Higgs vev). In the same way, M is the "bare" weak vector boson mass. 
The parameter f3 is such that (0|/i|0) = to all orders in perturbation theory. Since we 
are going to neglet the loop corrections, we set P = hereafter. Note also that the poten- 
tial in the Lagrangian (12. ip yields a partial spontaneous symmetry breaking (SSB). In fact, 
ground-state conflgurations are only invariant under U{l)em C SU{2)l ® U{1)y- 

The starting point of our work is the assumption that the Higgs fleld is a time- dependent 
"classical" homogeneous condensate, v = v{t). Working in the unitary gauge (^° = 0), the 
Lagrangian (12. ip becomes 

Ch = \v' + - \' + (2.4) 

where a dot denotes the derivative with respect to t and Cv,h is a polynomial in h and its 
derivatives. We then introduce a dimensionless real scalar fleld $ such that 

M 

v = 2 — ^, (2.5) 
9 

and 



V2 



<l>(t) + h{x) 



(2.6) 



9 

in order to study the fluctuations around the constant value <l>^ = 1, which is at the basis 
of the SSB mechanism. On neglecting all terms involving h (which will be analysed in 
Section [71), the Euler-Lagrange equation for $ becomes 

/\f2 f,,f2 

$_i!l^i$ + i!<fi$3^Q^ ^2.7) 
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Figure 1: The Higgs vev potential V{^) = $^ — 2$^. Straight horizontal lines correspond to 
c — —0.9, —0.5, 0.5, 0, 1, with c proportional to the Higgs vacuum total energy. 

which can be put in dimensionless form by rescaling the time as r = MHt/2 (primes will 
denote derivatives with respect to r), 

$" - 2<I>(1 - = 0. (2.8) 

After a simple integration we obtain 

$'2 _ 2$2 + $4 = $'2 + = c, (2.9) 

where the integration constant c e M is a function of the initial conditions. From a physical 
point of view, c is proportional to the total vacuum energ 30 and V{^) is the potential in 
which this vacuum is "moving" (see fig. [1]) . 

The equation (12.91) can be formally solved by separating variables, which yields 



V -'^'^ + 2<i>2 + c y$o 

Since we are interested in real solutions of Eq. (12. 9p . we must restrict to the case 

(_$4 + + c) > 0, (2.11) 

which corresponds to positive kinetic energy for $ (see eq. (12. 9p ). For any c, the stationary 
points of /($) are given by $ = 0, ±1, for which the potential V{^) has a maximum and 
two minima rispectively. Moreover, note that for c = and c = — 1 the integral in Eq. (12.101) 
is not defined at these points. We shall therefore consider the following cases: 



^It is easy to see that the integral over space of the Hamihonian density obtained from (|2.4p gives the 
total Higgs vacuum energy £■$ = Mhc/X. We shall thus refer to c as the "Higgs vacuum total energy" for 
brevity. 
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1. c < —1. Since c = — 1 corresponds to the absolute minimum in the vacuum total 
energy, this case is not physically relevant; 

2. — 1 < c < 0. The system exhibits closed trajectories in phase space and a periodic 
evolution; 

3. c = 0. The motion is along the separatrix and the system leaves the equilibrium point 
$ = with an exponentially growing velocity; 

4. c > 0. The system again exhibits closed trajectories in phase space but with periods 
longer than those in case [21 

It can further be shown that the solutions for all four cases can be connected by analytic 
continuatioiil. Therefore, in what follows, we shall use the simplest form of such solutions 
[20] (which corresponds to c > 0) for the sake of simplicity, namely 

<l>(r) = L^ cn(i?„r,m), (2.12) 



where Lm = J ^^"^ = J 2m-i ■ "^^^ parameter m is given by 



i + VTT~c 

m = , , (2.13) 

and is a function of the initial conditions because so is the vacuum energy c. The function 
cn{z,m) is called Jacobi Elliptic function and has the following properties |21j : 

periods: 4K, 2K + 2iK' (2.14) 
zeros: (2/ + 1)K + 2mK' (2.15) 
poles: = 2ZK + {2n + l)iK' , (2.16) 

where /, n are integers and K(m) is a special case of the complete elliptic integral of the first 
kind F{9, m), 

K(m) = F (^|, K'(m) = F (^^, m') , (2.17) 

with m' = 1 — m. From f l2.14p we deduce that the period along the real (imaginary) axis of 
the function in Eq. fl2.12p is given by 

{2v/2V2m- lF(f,m) , for m < 1 (m > 1) 

(2.18) 
v^V2m - 1 [F (f , m) + (|, 1 - m)] , for m > 1 (m < 1) . 

It is also important to note that the solution (I2.12p describes the evolution in time of the 
Higgs vacuum $ for every initial conditions ($0) '^'o) only after a suitable time shift (so that 
$'(r = 0) = 0). 



^ These solutions are well known in the literature for the theory in (1 + 1) and (3 + l)-dimensions 
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3 Fermions in a time dependent Higgs vev 



In the Standard Model the couphng of a generic fermion field ip to the Higgs scalar field H 
and gauge fields is described by the Lagrangian density terms [H] 

+^i^Rr (d^ - ^a'^B^j i^R-G {iJLHi^R + ^rH^iIjl) , (3.1) 

where ipi is a left handed isospin doublet, ipR ^ right handed isospin singlet and G the 
coupling constant between the Higgs boson and the fermions. Neglecting the gauge fields, 
in the unitary gauge (12. 6p . the previous Lagrangian becomes 



G 2M - G - 



with tjj = ^pL + il>R. Introducing the "mass parameter"|£l 

mfit) 



G 2M^, , G Mh^, , 



V2 g V2V\ 

with given by Eq. fl2.12p . the equation of motion for ip is given by 



G 

i-f^df, - mf{t) --^h{x) 



ij{x) = 0. 



(3.2) 



(3.3) 



(3.4) 



For our purposes, the last term above can be neglected with respect to the second term 
proportional to the Higgs condensate and one finally obtains 



[i-f>'d^-mfit)]ij{x,t) = 



(3.5) 



which is a Dirac equation with a time dependent mass. In a similar manner, one finds that 
the equation of motion for ip = ?/^^"7° is the Hermitian conjugate of (13. 5p . 

The spinor field is normalized in such a way that J d^xip\x,t)ip{x.,t) = 1 and, in 
the Heisenberg picture, it becomes a field operator with the usual anticommutation rules 

{ip{x,t),ip'^{x' ,t)} = 6{x — x'). We can then expand ip as 



ip{x) 



(27r)3/2 



5^e*^" \Us{k,t)as{k) + Vs{-k,t)h\{-k) 



where s = ± is the helicity, 

{as{k),al,{k')} = {h,{k),hl,{k')} = 6ss'S{k-k') 



(3.6) 



(3.7) 



Since there are no stationary states in a time-dependent external field <l>(i), the mass is strictly speaking 
ill-defined. We shall however refer to the function mf{t) as a time dependent mass. 
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and 



C/t(fc,t)f/,,(fc,t) = V,Hk,t)VAk,t) = 5s 





(3i 



ulik,t)vAk,t) = v,^ik,t)uAk,t) 

The vacuum state |0) is as usual defined by the relations 

as{k)\0) = bsik)\0) = 0. 



(3.9) 



With no time dependence in the theory (i.e., for mf(t) = nif constant), Us(k) and Vs{k) 
would be eigenstates of the operator 7 ■ k with eigenvalues rrif and —rrif respectively. The 
spinors in momentum space Us{k,t) and Vs{k,t) satisfy the charge coniugation relation 



^Uj{k,t) = Vs{-k,t), 

with ^ = z7o7J^. 

It is now convenient to introduce two new scalars defined by 

Usik,t) = [^7% + 7-^ + my(t)] X^+)(t)M, 
Vs{k,t) = \i-f^do-^-k + mf{t)] Xi-\t)vs, 



(3.10) 



(3.11) 



where 



Us 



(3.12) 



with xlXs = 1 and rjs = —ia2X-s, are eigenvectors of 7° with eigenvalues +1 and -1, re- 
spectively. Note that = v^g and Eq. (I3.10p is thus identically satisfied. With these 
notations, Eq. 03.51) yields 

W + Hit) T ^mf{t)] Xf )(t) = 0, (3.13) 
which is of the harmonic oscillator type with the complex and (doubly-)periodic frequency 

nl{t) T irrifit) = P + m){t) T irrifit). (3.14) 

Let us now assume that the Higgs vev remains constant and equal to $(0) for t < 0. 
Consequently, the mass mf(t) will also be constant at negative times and one just has plane 



*We are using the gamma matrices 



I 








-I 














where it-', j — 1, 2, 3 are Pauh matrices and I is the 2x2 identity matrix. 
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waves for t < 0. The evolution for t > is then obtained by imposing the following initial 
conditions at t = (jfl 



Xf)(0) = {21]fc(0) [nk{0) + mf{0)]} 



'1/2 



(±), 



(3.15) 



These together with Eq. (13.131) give 

At>(<) = (A-r'(*))' 



(3.16) 



SO that positive and negative energy modes are not independent and we shall then consider 
mostly the equation for xj^~^^ for simplicity. 



It can be showed that if fi{t) and f2{t) are two arbitrary solutions of Eq. (13.131) with the 
sign (+) (or, equivalently, with (-)), the quantity 



(3.17) 



is a constant of motion and one can then prove the stability of any arbitrary solutions [22]. 
Finally, note that if fi{t) = f2{t) = xl~^\t) the relation (I3.17P takes the form 



(+) 



+ .m,(t) (xw*xw-xr)*xw) = i. 



(3.18) 



which is also a consequence of the fact that Us{k, t) and Vs{k, t) are evolved by the Hermitian 
operators i'j'^do ^ ^ ■ k — mf{t). 

3.1 Fermion solutions and physical quantities 

The Hamiltonian operator for a fermion field can in general be written as 



n{t) = i / d^xij\x,t)'ip{x,t). 



(3.19) 



Inserting the expansion (13. 6p and using Eq. (13. 8p . this becomes 



n{t) = d^kJ2{ [iUlik, t)U,{k, t)J a\{k)as{k) + ^/(-A?, t)K(-A?, t)J hs{-k)h\{-k) 

+ \%Vl{-k,t)tj,{$A h,{-k)a,{k) + \iUl{k,t)Vs{-k,t)\ al{k)hl{-k)] , (3.20) 



■''We have set the momentum k = (0, 0, k). 
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where we have integrated on x and one of the momenta. Taking into account Eq. fl3.1ip . we 
end up with 



n{t) = d^kJ2^k{t)\E{k,t)\al{k)a,{k)~bs{-k)bl{~k) 



s=± 



where 



+F{k, t)bs{-k)as{k) + F*{k, t)al{k)bl{-k) 



F{k,t) 



au' 



(3.21) 



(3.22) 



(3.23) 



E'{k,t) + \F{k,t)\' 

From Eq. f l3.22p . using f l3.15p . it is possible to see that E{k, 0) = 1 and F{k, 0) = 0, therefore 
Tiit = 0) is diagonal. In fact, we have assumed that for t < the vacuum is constant and 
there is no explicit time dependence in the theory. 

The Hamiltonian 03.211) can be diagonalized at every time using a canonical Bogoliubov 
trasformation ^23j. As a matter of fact, the necessary conditions for this kind of diagonal- 
ization are ensured by the relations [T3] 



(3.24) 



iV}{-k,t)Vs{-k,t) = iUl{k,t)U,{k,t) = Qk{t)E{k,t) 

iUl(k, t)Vs{-k, t)j * = tV,^{-k, t)Us{k, t) = nk{t)F{k, t). 
We now introduce time- dependent creation and annihilation operators, 

as{k) 



as{k,t) 
blik,t) 



a{k,t) (3{k,t) 
-(3*{k,t) a*{k,t) 



as{k) 
bli-k) 



A{k,t) 



bli-k) 



(3.25) 



and the condition that this be a non singular canonical trasformation requires that ^ is a 
special unitary matrix, 

laik.t)]'^ + \(3{k,t)f = 1. (3.26) 

We have thus shown that SU{2) is the dynamical symmetry group for fermion creation in a 
homogeneus non stationary scalar field and the time-dependent vacuum |0)( (see below) is 
a generalized coherent state built on this group [21] . 
The hamiltonian also takes the diagonal form 



nit)= / d^kJ2^k{t) \dlik,t)dsik,t) -bs{-k,t)bl{~k,t) 



s=± 



(3.27) 



^This is a consequence of p.lSp . 
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if the coefficients of the canonical trasformation are such that 

(3.28) 

a{k,t) _ F{k,t) _ 1 + E{k,t) 
(3{k,t) ~ 1-E{k,t) ~ F*{k,t) ' 

which are indeed compatible with the condition fl3.26p thanks to fl3.23p . 

It is now possible to use the operators a e 6 to define time dependent Fock spaces, each 
of them built from the zero (quasi)particle state at the time t, ds{k,t)\Ot) = = bs{k,t)\Ot), 
which, at t = 0, are equal to as{k)\0) = = 6^(^)10). These relations mean that a quantized 
fermion field interacting with a classical external field $(t) can be represented at every time 
as a free field, with a corresponding redefinition af the particle concept and vacuum state. 
Moreover, one can show that diagonalizing the Hamiltonian f l3.2ip is equivalent to finding 
exact solutions of the Heisenberg equations of motion and all the matrix elements (expecta- 
tion values of physical observables) of interest can be written in terms of the coefficients of 
the Bogoliubov trasformation (13.251) [13]. 

For example, the vacuum expectation value of the (quasi-)particle number operator is 
given by 

Nkit) = (0|at(fc,t)a,(fc,t)|0) = \P{k,t)f {OMMmO) = \P{k,t)f 5{6), (3.29) 

where we have used Eqs. (13.250 and (13. 9p . From the previous relation we see that the 
number of created (quasi)particle pairs is spin-independent because the homogeneus field $ 
is isotropic [25]. If we put the system in a finite volume V, we must replace 6{0) in (I3.29P 
with = 1. The (quasi)particle density at the time t is thus given hyij 



= (oIt^ > : / (^s^'^wio) = (^/ d'kiomm = i / dkem,t)\\ (3.30) 



which is different from zero whenever the Hamiltonian is not diagonal in terms of the oper- 
ators a and b. The occupation number of fermions created with a given momentum k will 
be Ukit) = \P{k,t)\'^, and the condition (13.260 ensures that the Pauli principle is respected 
at every time [iHl EE] . 

In order to implement numerical methods, it is useful to cast some of the previous ex- 
pressions in dimensionless form. We thus introduce the following quantities: 



-Mh^ _ 2k _^G^ m) 



where m/ = m/(r — oo) = mf{Q). On further multiplying by A/ Mfj, Eq. (I3.13P takes the 
dimensionless form 

Xf^"(r) + [fi:2 + g$2(r,m)TZv^$'(r,m)]xfV) =0. (3.32) 



^Of course, the same result holds for antifcrmions. 
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Moreover, if we define the dimensionless frequency 

= ^ = V'^^^Ti^ , (3.33) 

the initial conditions (13.151) become 

(0) = { 2 ^,(0) [^,(0) + $(0, m)] Y^'^ 

(3.34) 

XF(0) = T^CU.(0)X(^)(0). 

Note that, despite this is not exphcitly indicated, the frequency and the initial conditions 
are also functions of the amplitude m of the elliptic function, that is of the vacuum energy 
c (see Eq. f l2A3D ). On using (I32HD, fl3:22D and (EUl) we finally obtain 

2 t^«;(r) 2u;«;(rj 

which gives the occupation number for every mode k as a function of the solutions = 
of Eq. (I3.32p . Analogously, the (dimensionless) energy density will be 



Pv(^) = j dKK ujf,{T)n^{T). (3.36) 

Note that n^^O) = p,/,(0) = thanks to the initial conditions fl3.34|) . 

A very useful result follows from the periodicity of the vacuum □, $(t) = $(r + T), which 
remarkably simplifies the evaluation of the occupation number and shows, although in an 
approximate way, its explicit time dependence [22]. If we define 

= uj^{t -oo) = uj^{0), (3.37) 

an approximate expression for nf^{T) is given by0 

, , {Im[xi'^(T)l}2 2 / r\ 2, , , , 

Mr) = ^ ■ 2,1 N sin^ = F.sin^ (u^r) , 3.38 

ui sm^(4) V TJ 

where Xk\t) satisfies Eq. (I3.32p with initial conditions xjc^\o) = 1, xj^''' (0) = and is 
such that cos{di^) = Re[xi^''(T)]. According to Eq. (13.381) the number density of fermions 
produced depends periodically on r for all k. So, on average, this density does not depend 
on the time during which the external field is turned on. The physical meaning of this result 
was pointed out by V. S. Popov through a quantum mechanics analogy [22]. Finally observe 
that in our case the time scale is fixed by the factor 2 / and if Mh ~ 10^ GeV then ["1 
2/Mh ~ 1.3- 10-26 s. 



^This definition is not restrictive since we only need either xi'^'^ or xi to calculate n,^{T). 
^We clearly refer to the period of $ along the real axis (see Eq. (|2.18p ). 

^°We shall see in Section [5] that Eq. (|3.38p is actually exact at t = nT for any positive integer n. 
"We recah that 1 GeV"! - 6.582 • IQ-^^ s for ?i = c = 1. 
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3.2 Band structure 



A non adiabatic quantum effect arises from tlie explicit time dependence of tlie frequency 
in Eq. f l3.32p . and tliis leads to tfie production of particles. When the time dependence is 
periodic, one usually speaks of parametric resonance. It is then clear that the quantity q has 
the role of a resonance parameter due to the fact that the time-dependent terms in Eq. (13.321) 
are proportional to q. 

We are interested in values of q and k which give solutions of the mode equation associated 
to particle production, identified by a mean occupation number (see Eq. (13.381) ) 

= {n^)r = (3.39) 

different from zero. The result is shown in Fig.[2]in which every peak corresponds to ra^ = 1/2: 
note the band structure in the plane (g, k^). The left plot in the upper part of Fig. [2] shows 
the first and second bands while the right plot displays bands from the second to the fourth. 
Moving along a band, oscillates between and 1/2. Moreover, the bands get narrower 
with increasing for a given value of q and after several bands they shrink to a negligible 
width, as we show in Fig. [S] with a plot of the mean occupation number as a function of 
for the bottom quark. 

To obtain the curves in the lower graphs of Fig. [2] we have studied the mean occupation 
number along straigh lines = dq, with d constant. First of all, for ^ El the peaks 
are located on the q axis at qn{'d ~ 0) = n^, with n a positive integer. We then made the 
ansatz qni^) = n^Vi^), with r(0) = 1, for the position of the peaks g„ along the generic 
line = dq. A numerical interpolation starting from the analysis of the first five peaks for 
different values of t? has then yielded 

m = (3-40) 
for which the parametric equation of the lines in Fig. [2] are given by 

g„(^9) =r(t9)n2 

(3.41) 

On replacing d from the first into the second equation, we finally find that the production 
rate must take its maximum values at points in the plane (g, k^) that satisfy the relation 

t^l + q^n\ (3.42) 

with n a positive integer. Note that the above relation represents a very good approximation 
in the regime of small oscillations around the static Standard Model Higgs. 

The band structure we just described will lead to a preferred production of non-relativistic 
particles, that is particles with small momentum compared to their mass. For example, for 

""^^Note that if i9 = 0, = and n^^, — with no production. 
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q 



4 6 8 10 12 14 16 18 

q 



Figure 2: Production chart for fermions coupled to the Higgs. Upper graphs show the 
occupation number as a function of the parameters q and . Lower plots show in details 
the lines in the plane (g, k^) along which the mean occupation number takes its maximum 
values. Note that for k ~ we have q^ — v?. 

fermions with a mass smaller than Mnjl^ the production is mostly driven by the first band 
in = 1) and Eq. fl3.42p yields a typical momentum fc^ ~ — ^"j which is smaller than 

rnj for the particles we consider in the following. On the other hand, the production of more 
massive particles will be caused by higher order bands (so that n"^ — q > 0) and is normally 
suppressed. 

4 Bosons in a time dependent Higgs vev 

The Lagrangian which describes the coupling between the Higgs field and the vector bosons 
is the same as in Eq.f l2.1l) . but with gauge covariant derivatives and kinetic terms for the 
gauge fields, 

Ch-b = D^H^D^H - f,'H^H - ^ {H^H)' - ^F;, - \g,. (4.1) 

13 
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Figure 3: Mean occupation number -Fk/2 as a function of for the bottom and Mh = 
115 GeV (solid line), 200 GeV (dotted line), 500 GeV (dashed line). Note that the peaks get 
narrower with increasing k^. Similar plots are obtained by varying the resonance parameter 
q at fixed as for decreasing q the bandwidth shrinks. 



where D,, = d,. 



- - ig% ■ - = d^Al - d,Al + ge'^^A^^A^ and G^, = d^B, - d,B^. 

For the electroweak vector bosons we thus have 



+$2 I M^W^^'W- + -yZ^Z^ 



9_ 

M 



2 ^ 

Ml 



+ 



$/i MW+^VT" + —Z^Z^ .(4.2 



Taking the vev with the condition (0|/;.|0) = and neglecting the hack-reaction of /i, the 
classical equation of motion for Z^ is given by 

d^&'Z^ + M|$2(t)Z^ = 0. (4.3) 

The fields satisfy the same equation with Mz replaced by M. We therefore conclude that 
any component ^ of the vector field Z^ and any component W of W^^ satisfy Klein- Gordon 
equations with mass M\{t) = M|$2(t) and M'^{t) = M'^^'^{t) which, after performing a 
spatial Fourier transform, take the form 



^fc+[fc2 + M|(t)] Ja = 0, 

n + [k'' + M\t)]n = o, 



(4.4) 
(4.5) 



where ^{t) = {2n)-^/'^ J d^k e'''-^ ^k{t) and W{t) = {271)"^/'^ J d^k e''' '^Wk{t). Analogously, 
from the Lagrangian (12.11) in the unitary gauge (12. 6p and keeping only terms up to second 
order, one obtains the following equation for the Fourier modes of the quantum Higgs field 



(4.6) 
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where 



2 ^ ' 2 



(4.7) 



We are interested in the case when the Higgs vev $(t) oscillates near one of its two 
absolute minima (see Fig. [T]). We therefore take the limit m — > +oo (c — —1) in Eq. 02.121) 
and obtain 

$(r) ~ 1 - — sinV + ( At I . (4.8) 

Moreover, from now on, we shall for simplicity denote with any bosonic mode, that is 
can be either #^ or h^. Eqs. (14. 4p . (14. 5 p and (14.61) then become 

>;"(r)+^^(r)r.(r) = 0, (4.9) 

with 

^l^r) = K' + q(l-^sin^T) , (4.10) 

where we have neglected terms beyond the first order in 1/m and introduced the dimension- 
less constants 

2k ^ Ml , , 

and ,.4-jj|, (4.11) 

with My = Mz, M or Mh and cry = 1 for the gauge fields and 3/2 for the Higgs. Note that 
Eq. (13.321) for the fermion modes takes the same form as Eq. (14.91) when the imaginary part 
of the fermion frequency can be neglected and that the period of ^^(r) is vr. 
If we define 

-ip-'(-a- ("^) 

and change the time to rj = 2t = Mnt, we obtain 

F;(r7) + (a + 2ecosr/)r«(r7) = . (4.13) 

When the total vacuum energy c is near the minimum value c = — 1 (for which $ is constant), 
W^'^, and h therefore satisfy a Mathieu equation. 

The Mathieu equation is a linear differential equation with time periodic coefficients and 
is described in the general framework of Floquet theory [2]- An important result from this 
theory is the existence of stable solutions only for particular values of a and e (see Fig. Hj). 
For our purposes, the relevant solutions which lead to an efficient particle production in the 
quantum theory are however those which show an exponential instability of the form [iiif* 
is known as Floquet index or growth factor) 

~ exp {ii^:^r^/2) (4.14) 

and appear within the set of resonance bands of width Ak*^"^ labelled by the integer index n. 
Using Eq. (14.120 one can then map the stability chart of Fig.lHinto the plane (g, k^) and finds 
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Figure 4: Stability regions for the Mathieu equation. Stable solutions have a and e in 
the unshaded regions. When e = 0, instability bands cross the a-axis at = with 
n = 0,1,2,.... 



that the instabilities also correspond to an exponential growth of the occupation number of 
quantum fluctuations, 72^ oc exp{jjl^\) (see Eq. fl4.18p below), which can be interpreted as 
strong particle production. Stable solutions also lead to particle production, although with 
no exponential growth in time, and thus resemble the fermion case [22]. 
If c ~ — 1 then m ^ 1 and e ^ in FigHl One then finds 4a ~ n^, or 



which, for a given q, gives the value of around which the n-th resonance peack is centred. 
Further, the physical constraint > implies q < 2n^ m/ (2m — 1), so if g < 1 {Mh ^ 2My) 
all the resonance bands contribute to production, otherwise at least the first band {n = 1) 
is not available. In particular, the first band is never available for the Higgs particle and 
its production can therefore be negligible with respect to that of gauge bosons if the Higgs 
mass is sufficiently large. 

4.1 Boson occupation number and energy density 

The same analytical steps of Section [3] allows one to relate relevant physical quantities to 
the solutions of the equation of motion also in the case of bosons. We again take initial 




n = 1,2,3, ... , 



(4.15) 
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conditions corresponding to positive frequency plane waves for r < 0, 



(4.16) 

The occupation numbers for bosons will then be given by [I3] 

n.{r) = -^Kn^\Y.\'-l, (4.17) 

where gj^ = ?a7^(0) and note that 11^(0) = thanks to eqs. (14.161) . 

Like for fermions, it is possible to find an analytic approximation for the boson occupation 
number [TH] . 

= sinh^(;U«;r) , (4.18) 
with the Floquet index, or growth factor, /i^ given by 

cosh(/x,T) = Re[y«(T)] , (4.19) 

in which Y^^^ is the solution of Eq. (14. 9 p with initial conditions Yk^\o) = 1, Fi^-* (0) = 
and T = vr is the period of Wi^{t) (the period of the Higgs vacuum). Using this approximate 
expression, the boson energy density can be written as 

PB(r) = ^ J dK k^w^{t) smh'^{n,,T) , (4.20) 

where the tilde on p is to remind that this is a dimensionless quantity (like the fermion 
analogue in Eq. (I3.36P ). 



5 Fermion production 

We shall now use the numerical solutions to Eq. (13.321) to evaluate the fermion occupation 
number in Eq. fl3.35p and compare it with its analytical approximation in Eq. (I3.38p . 
In order to illustrate the magnitude of the effect, we shall consider three possible values 
for the Higgs boson mass which result from different experimental lower or upper bounds, 
namely Mh = 115 GeV, 200 GeV and 500 GeV [H]. In this Section, we shall not take into 
account the back-reaction of the produced fermions which is treated later. 

Let us begin with Fig. El which shows the time evolution of the occupation number 
and its enveloping for the top quark in the case Mh = 500 GeV and two different values 
of the momentum k, and Fig. El which shows the same quantities for the bottom quark in 
the case Mh = 200 GeV: 
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Figure 5: Time evolution of the occupation number and its envelope for the top 
quark and Mh = 500 GeV. Left graph: = 1 — g ~ 0.51 on the first resonance band, 
Co = — 1 + lO"'^ (solid line) and c = — 1 + 4 ■ lO^'^ (dashed line). Note that and exactly 
coincide and their period scales as (c + 1)^/^. Right graph: = g/10 ~ 4.8 ■ 10"^ is outside 
resonance bands and (solid line) differs from (dashed line). 

a) For K, on the (first) resonance band both functions and for the top reach the 

maximum allowed by the Pauli blocking and are remarkably indistinguishable (left 
graph in Fig. [5]). In fact, their difference remains of the order of 10"* for both values 
of the initial background energy c = c(r = 0) and we do not show it. For the bottom 
{i.e., a smaller value of q at the resonance with respect to the top's) instead, the two 
functions differ slightly and always remain smaller than one (see Fig. [6]). In both cases 
however, from the numerical simulations one can infer a scaling law for the period of 
the occupation number with respect to the initial vacuum energy c. If Cq is a reference 
energy, one has 

Uc)=Uco)(^^y , (5.1) 

which holds for all values of k on the n-th resonance band provided c is small enough 
that only small oscillation of the background $ are relevant (see case [2] in Section [2]). 
Note also that ^ T for the top (we recall that the Higgs vev period is 2T/Mh ~ 

10-26 s). 

b) For values of k not on a resonance band (right graph in Fig. E]), the production is of 

course much damped and the exact occupation number shows relatively high frequency 
oscillations, with period comparable to the Higgs few's, which are modulated by the 
function with a period usually larger than T. It is therefore clear that the time 
average over these relatively high frequency oscillations can still be related to the 
enveloping as 

1 

n^ir) = {n^{T))r = f c?C^«(C) ^ n^{r) , (5.2) 
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Figure 6: Time evolution of the occupation number for the bottom quark, Mh = 200 GeV, 
= 1 — g ~ 0.9982 on the first resonance band and Cq = — 1 + 10"'^. Left graph: (dashed 
hne) and its envelope (solid line). Right graph: difference between and n^. 



where is given in Eq. fl3.38p . 

We emphasize that the fermion occupation number is always periodic in time and the function 
in Eq. (13.381) exactly equals rif^^pT), with p a positive integer, regardless of the fact that k 
and q are on a resonance band or not (see Fig. E]). This is a general behaviour which still 
holds for very small values of q, for example q ~ 10~^° or g ~ 10~^^, corresponding to the 
scales of the electron and neutrino mass (< 1 eV) respectively. These cases are however very 
difficult to treat numerically because one should integrate for very long times (the period of 
the function in Eq. (13.381) grows with the inverse of the fermion mass). For this reason we 
have chosen to show the occupation number only for quarks like the top and bottom. 

We are now interested in the production probability for a given fermion, i.e., in the 
occupation number summed over all k's for a given value of q (see Eq. (13.311) ). The 
generally oscillating behaviour of leads us to estimate this probability in time as the 
mean number (15.21) summed over all the momenta k, 



4 7r3 / 2 



(5.3) 



As we see from Fig. |3l the mean occupation number as a function of k at fixed q is 
significantly different from zero only near the values = k\ around which the production 
peaks are centred. We can therefore estimate 



rip 

4 7]-3 

n=l. 



^ rip 
2^^ 



(5.4) 



Supp{Pn) 



n=l 



where Supp{Pn) is the interval around in which is significantly large. The effective 
width of the peak Ak^ on the n-th band is determined by estimating each integral in the 
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Mh = 115 GeV 


Mh = 200 GeV 


Mh = 500 GeV 






Q 






Q 






Q 


t 


< 10-10 


9.26 


t 


1.6 - 10-5 


3.06 


t 


4.3 - 10-^ 


0.48 


b 


8.7-10-5 


5.46 - 10-3 


b 


5.1 - 10-5 


1.80 - 10-3 


b 


2.4 - 10-5 


2.90 - 10"^ 


d 


1.6 - 10-^ 


1.38 - 10-^ 


d 


8.2 - 10-^ 


4.54 - 10-9 


d 


3.3 - 10-s 


7.28 - 10-1° 


T 


3.9-10-5 


9.54-10-^ 


T 


2.4-10-5 


3.16 - 10-^ 


T 


8.8 - 10-6 


5.06 - 10-5 



Table 1: Production probability of various quarks and of the lepton tau for three different 
values of the Higgs mass and Higgs total vacuum energy cq = — 1 + 10"^. 



above expression numerically. Further, we found that 

A^lic) A^lico) (^)' , (5-5) 

a scaling behaviour also shared by the bosons (see Section [6]). 

The occupation number in Eq. (15. 4p for the top, bottom, down and for the tau are shown 
in Table [T] for three different values of the Higgs mass and Higgs total vacuum energy 
Co = — 1 + 10-3. Note that the production is mostly generated from the first band and, 
if g < 1, the production probability is larger than for g > 1 since for the latter case only 
the bands of order greater than one can contribute. So, recalling that q = 4m'j/Mfj, 
we conclude that in the Standard Model all the fermions except the top have a greater 
production probability for a "light" Higgs. Moreover, within a given band, the production 
is still strongly affected by the resonance parameter q. On varying q along a given band one 
observes a change in the value of the momentum k of the produced particles (see Eq. (I3.42p ). 
in the width of the band and, consequently, in the number density of produced particles (see 
Figs. [2] and [3]). These facts lead to the values of (as defined in Eq. (15. 4p above) presented 
in Table [H Note finally that the bottom is the "dominant" (with the highest production 
probability) fermion for Mh = 115 GeV and Mh = 200 GeV, while for Mh = 500 GeV the 
top production is more probable. 

We remark that the average occupation number also scales with the energy according 
to the relation 

^^y~(l + c)t, |c+l|<10-3, (5.6) 

where n is the order of the "dominant" peak (the one which contributes most to the produc- 
tion). So if one is interested in the production probabilities at an arbitrary but always very 
small energy, the values in Table [H calculated for cq = — 1 + 10"3^ must be multiplied by the 
factor [(1 + c)/(l + Co)] 2 if the production is associated with the first band. For example, 
this is true for the top if Mh > 2 mtop and all other fermions for any Higgs mass. Instead, 
in the region Mh < 2mtop (of phenomenological interest), the top may be produced only 
starting from the second band and the relevant scale factor is (1 + c)/(l + Co). 
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Finally, we evaluate the fermion energy density using Eq. 03.361) . The integral over k can 
be calculated in the same way as in Eq. fl5.4p and including just the dominant peak for the 
production. This approximation is better the lower the Higgs vacuum energy, since the peak 
amplitude decreases proportionally to this energy, we thus find 

1 /" F If 

P^{t) = (i^K t^«,(r)-^ sin^ (z/^r) ~ — / c?^/tt^«:(r) n« sin^(z/«;r) . (5.7) 

Supp{Pao,n) 

Around the dominant peaks given by Eq. fl3.42p we have 

^ 1 , (5.8) 
and we can therefore write the energy density fl5.7p as 



^ , , sin (r) /" ,o , , sin (r) vA . 9 / \ ,^ ^^ 

P^i'^)-^r^ d^t^uj^{T)n^^ l^^^A/s:>«„(r)n,„ , (5.9) 

which significantly simplifies the evaluation of the integral over k. One can also find a scaling 
law for the energy density so obtained, namely 

p^~(l + c)?, (5.10) 

with the same prescriptions as for n^. 

From Eq (12.41) . the vacuum energy in a volume 8/Mfj can be calculated and put in the 
usual dimensionless form as 

£,^^£, = lc. (5.11) 
The Higgs vacuum energy density available for the production is thus given by 

=S^- 4=1 = Y (c + 1) • (5.12) 
A 

In the Standard Model A appears as an arbitrary parameter and there are presently no 
experimental constraints on it. There are however good reasons to believe that A < 1 but 
we shall consider in all our calculations simply the value A = 1. Taking into account the 
dependence of (p^) and A£^$ on c, the fraction of total energy absorbed by a given fermion 
scales as 

R(c) = ^ = Rico, X) f^y^', i?(co,A) = ^(l + Co)-"/^(p4co)) (5.13) 

An example of this quantity is plotted in Fig. [71 Notably, the maximum value is reached 
before the end of the first vacuum oscillation. 

Using the results shown in the Table [Tj for each included fermion we have selected the 
values of the Higgs mass which give the greatest production probability and the corresponding 
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Figure 7: Fraction of Higgs vacuum energy R absorbed by the produced top for = 
500 GeV and cq = -1 + lO'l 





R 




Mh 


t 


1.1 • 10-1 


0.48 


500 GeV 


b 


2.2 • 10-2 


5.46 ■ 10-3 


115 GeV 


d 


4.0 ■ 10-^ 


1.38 ■ 10-^ 


115 GeV 


T 


1.0 ■ 10-2 


9.54-10-4 


115 GeV 



Table 2: Fraction of Higgs vacuum energy absorbed by various quarks and the tau lepton 
for the particular Higgs mass (out of the three considered in Table [1] with cq = —1 + 10-^) 
which gives the highest production rate. 

energy densities are shown in Table [2] in units of the Higgs vacuum energy density avaible 
for the production A£^$. Of course, the physical condition R < 1 limits the range of validity 
of this first approximation where production is not taken to affect the vev oscillation (no 
back-reaction), breaking energy conservation. For example, the value of R given in Table [2] 
for the top for cq = — 1 + 10-^ can be rescaled down to a minimum energy c ~ — 1 + 4 • 10-^, 
for our choice of A. 

When their energy density becomes comparable with the background energy, the pro- 
duced fermions are expected to back-react on the Higgs vev, thus affecting its evolution and 
eventually suppressing the parametric production of particles. We shall see this in detail 
in Section [71 For now we just note that the different scaling law for the top quark energy 
density as a function of c (when the Higgs mass is large) implies that back-reaction effects 
will become important sooner and one should therefore be aware that in this case, for very 
small Higgs energy oscillations, the top energy density produced will be strongly affected by 
these effects. 
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Figure 8: Floquet index for the first resonance band function of the Hide's mass for 

(dashed hne) and (sohd hne). In both cases /il^'* is maximum if Mh ~ 2 My and 
decreases for larger values. A factor vT+~c on the vertical axis is used in consideration of 
Eq. dSH) with Co = -1 + lO"^. 

6 Boson production 

In Section H] we reviewed the fact that the equation of motion for a bosonic field admits 
unstable solutions only for some values of the parameters q and k. The results of the 
numeric integration of Eq. (14.91) will be used to evaluate the occupation number for every 
mode Kn given in Eq. (14.151) . For both Zq and we have seen that, for a given value of 
Mh, the first band (n = 1) contributes to the production in an extremely dominant way 
with respect to the others since the Floquet numbers scale with the Higgs vacuum energy 
according to 



where we usually take cq = — 1 + 10^'^. Therefore, in the following we shall just refer to the 
first band and then do not describe the Higgs production (see Section H]), which could be 
trivially included. Moreover, motivated by the discussion at the end of Section HJ we shall 
first study the case Mh — 2 M^o = 182 GeV as the plot in Fig. [H] shows that this yields a 
more efficient production. 

Let us note that for the case of production dominated by the second band {i.e., for bosons 
with mass less than twice the Higgs mass, including the Higgs itself), the relevant Floquet 
index is much smaller. In particular, a plot similar to that in Fig. [8] would display a curve 
for \ which scales as (1 + c) instead of (1 + c)^/^, with maximum value on the vertical 
axis of the order of 10~^. 

Fig. shows the time dependence of the occupation number (I4.17P and its enveloping 




(6.1) 
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Figure 9: Time evolution of the number of produced with ~ 0.17 for two time 

ranges, Mh = 200 GeV and cq = — 1 + 10~^. Note that (sohd hne) matches (dashed 
hne) at the end of each background oscillation. 



approximation fl4.18p for the boson with a momentum on the first resonance band, = 
n\ ~ 0.17, Mh = 200 GeV and Cq = — 1 + 10~^. The exact n^, whose mean value grows 
exponentially with time, oscillates and coincides with at the end of each background 
oscillation. 

To evaluate the production probability {i.e., the total occupation number) for a given 
boson we must integrate the mode occupation number over a finite volume in momentum 
space (at fixed q). This very difficult integration is greatly simplified by the fact that the 
production only arises around the peaks in the (g, k"^) plane, that is for = given in 
Eq. (14.151) . Moreover, the production from the first instability band is the most relevant, 
so that the integration can be consistently restricted to an interval around ni (denoted as 
Supp{Pi), like for fermions), 

risir) = j^-^ j d^Kh^{T) ^ jdn'^ v^sinh^ (/x^r) . (6.2) 

Supp{P-,) 

We shall then consider two different time scales. For short times (/i^T ^ 1), the above 
expression can be approximated as 

UBir) ^ jdK^ Vk^ {fx^rf , (6.3) 

Supp(Pi) 

whereas for long times (/U^t ^ 1) we shall use 

nnir) -i^Jdn^^^ e^^'^" . (6.4) 

Supp{Pi) 
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Figure 10: Time evolution of the total number of produced (left graph) and fraction of the 
Higgs vacuum energy absorbed by (right graph) for Mh = 500 GeV and cq = — 1 + 10~^ (in 
the right plot, fractions larger than unity signal that the back- reaction cannot be neglected). 



The Floquet index as a function of has a parabolic shape around ki, 



(6.5) 



where AkI denotes the width of the first peak (the same notation we used for the fermions) 
which scales with the Higgs vacuum total energy exactly according to the same law (15.51) for 
fermions provided c < — 1 + 10^'^. 

For /i^T ^ 1, the integral in Eq. (16.30 can now be easily estimated and yields 



2 7r2 



(6.6) 



where Nb = Nb{ki, Anf) is a rather cumbersome expression which we do not show explicitly 
since it will not be used in the present paper. For long times, one can likewise estimate the 
integral in Eq. (16 ■4p using a saddle point approximation and obtain 



nB(r 



16 TT^ 



TT 



2t 



(6.7) 



The left plot in Fig. [10] shows the time dependence of the total number of produced 
Z^ following from (16. 7p for Mh = 500 GeV and Cq = — 1 + 10~^. This quantity grows 
exponentially in time and when the fraction Efl of Higgs vacuum energy converted into Z^ 
approaches unity (see right plot), the back-reaction is expected to become relevant. Indeed, 
values of that fraction larger than unity do not make sense and actually signal the complete 
failure of the approximation which neglects the back-reaction. Note that this occurs after 



-"^•^This fraction is defined in analogy with that for the fermions given in Eq. (|5.13p . 
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about a thousand background oscillations. Recalling the results for the top from Section El 
it is clear that, during the first few hundreds of vacuum oscillations the energy transferred 
to the top is larger than that absorbed by the Z^, and it seems that the back- reaction of 
the latter particles can be neglected at an early stage. The production of and their 
interaction with the time dependent Higgs background take the lead later and, with a good 
approximation, remain the only processes with significant effects (we shall have more to say 
on this in the next Section). 

7 Fermion and boson back-reaction 

The results presented so far have been obtained by neglecting the back-reaction of the pro- 
duced fermions and bosons on the evolution of the Higgs vacuum. This is a common approx- 
imation, since a more complete treatment of these effects is a very difficult task. We shall 
therefore begin by considering the simple case of one kind of fermion simply denoted by ip 
(which we shall identify with the top quark later) and one kind of boson (the Z^). 

The Lagrangian density for our system is given by (see Eqs. (12.41) . (13.21) and (14.21) with 
V = Mh^/V\) 

Cbr = '^h-^ V^^i^ $ - ^ g zf^z^ $2 , (7.1) 

where we have redefined the product (ipip) to make it dimensionless. The back-reaction effects 
can be studied to a good degree of accuracy in the Hartree approximation (see Refs. [151 E] 
for a discussion on this subject). Let us then consider the vacuum expectation value of the 
Euler-Lagrange equation for $ which reads 

$" - 2$(1 - $') + ^ v^(^^) - ^ q {Z^Z,)^ = . (7.2) 

The product (ipip) can be rewritten in terms of an integral in momentum space of the 
Bogoliubov coefficients (previously introduced to diagonalize the Hamiltonian operator for 
a generic fermion field) using the expansion (13. 6p . In this way one encounters ultraviolet 
divergencies and, in order to obtain a finite result, the operator ipip must be normal ordered 
to subtract vacuum quantum fluctuations [m [151 [29] . Moreover, we shall neglect any other 
renormalization related, for example, to perturbative quantum corrections. Note that after a 
Bogoliubov transformation the vacuum state \0r) is time-dependent and we have a different 
renormalization at every time. We thus define the (time-dependent) normal ordering of a 
generic operator as 

^r{O) = O~{0r\O\0r), (7.3) 

so that the vev of ipip will be given by 

{^P^) ^ {0\yy;(i,i,)\0) = (0|^^|0) - (0|((0.|^^|0.)) |0) 

= {0\^m - (0.|#|0.) . (7.4) 
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Using results from Section [HI we find 
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2 7r3 
-2n^(r) + 



X, 



2^3 



T - 



Im [X^{t)X':{t 



(7.5) 



As for the Z^, analogous prescriptions to those used for the fermion lead to 



d' 



K 



(2 



\Y r - 



-^^2 + g $2 



(7.6) 



Note that the two back-reaction terms (17. 5p and (17.61) vanish at r = by virtue of the initial 
conditions (13.341) and (I4.16p . as expected. The system of back-reaction equations is thus 



$" _ 2$(1 



XI + + g$2 - X, = 

I y; + (/s:2 + g<|.2)y; = o. 



(7.7) 



The last two terms in the first equation (Higgs-fermion and Higgs-vector coupling terms) 
depend on A, that is the strengh of the quartic term in the Higgs vev potential. In the 
Standard Model the coupling A is not fixed but it is possible to take < A < 1 if one 
considers the Higgs self-interaction as described by a perturbative theory. 

For the above system, it is easy to show that the (renormalized) total energy 



A 



f2 7r)3 +2 



d' 



K 



3 ^K. II'K 



is exactly conserved by virtue of the equations of motion themselves. In the above 

c(r) = [^'{r)f -2^\T) + ^\r) , 



(7i 



(7.9) 



the fermion (boson) occupation number tIk {tlk ) and frequency {zUf^) are given in 
Eq. fl335|) [Eq. KT7\l ] and Eq. (K33\l [Eq. fITTOD ]. respectively. 

The approximate evaluation of the integrals (17. 5p and (17.60 by restricting their integrand 
on the dominant band (see also [T^) is based on the results of Sections and El for fermions 
and bosons respectively, and proceeds in a similar way. We therefore assume the validity of 
scaling laws of the form given, for example, in Eqs. (15. 5p or (15. 6p and Eqs. (16. 6p and (16. 7p . 
These approximations are not valid for asymptotically long time evolution, when the dissi- 
pative dynamics takes over and very little particle production may still take place. In such 
a case, one may need different computational schemes, and probably a full lattice approach. 
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useful also to describe rescattering phenomena and estimate a possible termalization phase. 
Nevertheless, our approximation scheme has the virtue of been simple enough to lead to 
a system of differential equations of finite order which helps to grasp some aspects of the 
back-reaction dynamics and gains more and more validity in the limit A -C 1. 

In order to make the qualitative features of the system (17.71) clearer, we shall first study 
separately its behaviour for short and long times. In particular, as we have seen in the 
two previous Sections, the production of fermions is expected to dominate during the first 
few hundreds of Higgs background oscillations and in that regime we shall neglect boson 
production by simply switching its coupling off and set Mh = 500 GeV in order to have 
the maximum top production. At larger times the production of bosons overcomes that 
of fermions (limited by Pauli blocking) and we shall then neglect the fermions and set 
Mh = 200 GeV so as to maximise the production. In the last Subsection, we shall finally 
consider the entire system with both the top and for = 500 GeV in order to show 
the effect of the fermion back-reaction on the later boson production. This value of the 
Higgs mass, although not very likely, is chosen as it is particularly convenient to study the 
interplay between fermion and boson production. 

7.1 Fermion back- react ion 

As we mentioned above, for relatively short times, we at first consider the evolution of the 
Higgs background and one fermion, namely the top quark with Mh = 500 GeV. 

The left graph in Fig. [TT] shows the time dependence of the corresponding occupation 
number for the top and Higgs energy densities up to a time r ^ 3.14 • 10^ (~ 10^ T, where 
T is the Higgs period without back-reaction) when the boson back-reaction is expected to 
take over. The Higgs energy oscillates periodically in time and of course takes its minimum 
values when the top's occupation number is maximum. The right graph displays how the 
Higgs energy changes in time for different initial values of c, that is for different total energy 
(all curves assume A = 1). Note that the maximum fraction of Higgs energy converted 
into the top increases for decreasing initial total energy, as one would expect from the Pauli 
blocking. Fig.[T2]then compares the top's occupation number with and without back-reaction 
for c(0) = — 1 + 10~^. It is clear that the back-reaction in general suppresses the number 
of produced fermions and this effect is more pronounced for smaller total energy. We also 
remark that the above numerical solution conserves the total energy with an accuracy better 
than 1 part in 10^. 

To conclude this section, we would like to spend a few words about the behaviour of the 
system with respect to the value of A. Further numerical analysis shows that the fraction 
of energy absorbed by the produced fermions is proportional to A. Moreover, the system of 
non-linear differential equations (17. 7p seems to produce a chaotic behaviour only for A > 10^. 
In fact for A ~ 10^ the highly non-linear term proportional to (ipip) becomes of the same 
order of magnitude as the other terms in the first equation of (17. 7p . But at this level the 
approximations used in our equations to evaluate the back-reaction terms should break down. 
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Figure 11: Fermion back-reaction with Mh = 500 GeV. Left plot: Time evolution of the top 
occupation number (dashed line) and Higgs vacuum energy (solid line) for A = 1 and initial 
c(0) = — 1 + 10~^. Right plot: Comparison of the Higgs total vacuum energy for initial values 
c(0) + l = 10~^ (dotted line), 10~^ (dashed line) and 10~^ (solid line). A normalization factor 
of [c(0) + 1]^^ is used for the Higgs energies. 
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Figure 12: Fermion back-reaction: comparison of top occupation number with (dashed line) 
and without (sohd line) back-reation for c(0) = — 1 + 10^^. 

7.2 Boson back-reaction 

From the results of the previous Section, we know that the energy density of created bosons is 
comparable to the initial Higgs energy density after about 10^ vacuum oscillations (that is, for 
r > 3.14- 10'^). The production of bosons then overcomes fermion production (limited by the 
Pauli principle) and their back-reaction becomes the most relevant phenomenon. We shall 
therefore neglect the fermion contribution here and just consider the and = 200 GeV 
or Mh = 500 GeV. 

Even with the above simplification, it is still very difficult to solve the system (17. 7p and 
we need to employ yet another approximation. Since the integrand in Eq. (17.61) is sharply 
peaked around the centres of resonance bands, we estimate that integral in the same way 
we used to obtain the total number of produced bosons in Eq. (16.71) . as we have already 
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anticipated in the general discussion of this Section. 
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Figure 13: Boson back-reaction with Mh = 200 GeV. Left graph: occupation number 
(dashed hne) and Higgs vacuum energy (dotted hne) for an initial value of c(0) + 1 = 10~^. 
Right graph: Time variation of the Higgs vacuum energy for initial total energy c(0) + 1 = 
10""^ (dotted line), 10~® (dashed line) and 10~^ (solid line). A normalization factor of 
[c(0) + for the Higgs energies is always used. 
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Figure 14: Boson back-reaction: comparison of occupation numbers with (dashed line) 
and without (solid line) back-reation for c(0) = —1 + 10~^. = 200 GeV in left plot and 
Mh = 500 GeV in right plot. 

We can now solve numerically the system (17.71) and Fig. [13] shows some relevant results. 
Starting from the left, it is possible to see that the Higgs energy is minumum when the rate 
of production is maximum. Further, from Fig. [TH one sees that the back-reaction strongly 
suppresses the boson production. The right plot in Fig. [12] shows the amount of vacuum 
energy dissipated by back-reaction effects as a function of the time for three different values 
of the total (dimensionless) energy (which of course coincides with the initial Higgs energy). 
An important feature of this plot is that we have a "regeneration" of the vacuum energy, 
which now oscillates in time, and the number of peaks in a given time interval is seen to 
increase (that is, the periods of oscillations become shorter) for increasing total energy. In all 
cases, the values at the peaks slowly decrease in time and the periods of oscillations stretch. 
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One can therefore conclude that the system will evolve towards a complete dissipation of 
the Higgs oscillations. This is due to the presence of a dumping term [analogous to 1/a/t in 
Eq. (16.71) ] that appears in the back-reaction term (17. 6p after integrating in k. 



7.3 Fermion and boson back-reaction 




Figure 15: Time evolution of the occupation numbers for the top (dashed line) and (solid 
line) for A = 1, Mh = 500 GeV and c(0) = — 1 + 10"^. Left plot: Partial hack-reaction from 
the approximations of Sections 17.11 and 17.21 Right plot: Full back-reaction horn Section [7[3l 

We finally consider the whole system with the top and described by Eq. (17. 7p for 
Mh = 500 GeV. Again, we choose this value of the Higgs mass in order to have both fermion 
and boson production from the first resonant band. This corresponds to a particularly 
complicated scenario with large production of both the top and Z^, whereas a more realistic 
smaller value of Mh would mostly lead to the production of lighter quarks with smaller 
effects on the Z^ at late times. 

The coupled dynamics is then described in Fig. [15] in which we show on the right the 
occupation numbers for the top and Z^. This results has to be compared to the results of 
back-reaction computed in the "factorized" approximation used in Sees. 7.1 and 7.2 which 
we report together for convenience on the left, both computed at Mh = 500 GeV. We note 
that the top production is not particularly affected, whereas the Z^ occupation number is 
suppressed in the sense that it increases significantly only at later times. The time needed 
for the production of Z° to overtake the top is about one order of magnitude larger than 
without fermion back-reaction. 

As we have already observed, we are able to study the evolution of the system in a win- 
dow of time far from the asymptotics (presumably with an upper bound in time inversely 
proportional to A). Beyond that, the adopted approximations break down and the asymp- 
totic time behaviour of the system should be studied taking into account a more accurate 
description of the mode production (eventually in a full lattice approach) together with the 
rescattering contribution of the produced particles. 
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8 Conclusions 



We have considered oscillating solutions for the Higgs vev in the context of the Standard 
Model of particle physics and studied some resulting dissipative effects. In the Standard 
Model, in fact, the masses of fundamental particles depend on that vev and its oscillations 
can be viewed as a time-dependent renormalization of the particles' masses which leads to 
the production of fermion and boson pairs by parametric resonance. 

In the first part of the present paper, the back-reaction of the produced pairs has been to- 
tally neglected. In this approximation, particle production by the oscillating Higgs appeared 
to be very efficient. For fermions, the Pauli blocking constrains their occupation numbers to 
oscillate in time (see Figs. E]) about mean values smaller than one. From the entries given 
in Table [H one can see that such mean values strongly depend on the resonance parameter 
q, that is the Higgs and fermion masses [see the definitions f l3.3ip ]. Moreover, the particular 
form of the governing Dirac equation f l3.13p yields a significant probability of producing only 
those fermions whose momenta lie on well-defined bands in the (g, K^)-plane [see Fig. [2] and 
Eq. fl3.42p ]. Considering the masses of Standard Model particles, this also implies a larger 
probability of producing fermions with non-relativistc momenta. 

As for the bosons, their production is not constrained by any fundamental principles 
and the governing Mathieu equation f l4.13p (in the small oscillation regime) leads to occu- 
pation numbers which grow exponentially in time. Analogously to the fermions, the boson 
production only occurs on narrow bands in the (g, fi;^)-plane [see the definitions (14. lip and 
Eq. glSD]. 

In our analysis, we have regarded the Higgs mass Mh as an adjustable parameter, for the 
experimental data only place a lower bound on it. Consequently, we have given our results 
for different possible values of Mh and the initial Higgs energy c(0) [see the definition (12. 9p 
and Eq. (I7.9p ]. For natural values of these parameters, we have found that a significant 
fraction of the initial vacuum energy can be transferred to the fermions before a complete 
background oscillation (whose period we denoted as T, see Fig. [7]). The Pauli blocking then 
prevents the fermion production from increasing further. The bosons, on the other hand, 
grow exponentially in time but absorb a significant fraction of the Higgs vacuum energy only 
after about a thousand oscillations (see Fig. [10]). This feature seems to suggest that it is 
possible to analyse the back-reaction of fermions and bosons separately. 

In particular, we have studied the system of one fermion (the top quark) and one boson 
(the Z°) coupled to the Higgs field $ in the Hartree approximation [see Eq. (17. 7p ]. For such 
a system, the total energy (17. 8p is conserved and, for < r < 10^ T, the back-reaction of 
the produced Z° has been neglected. To make the problem more tractable, and avoid a 
complicated study of the system on a lattice of Fourier modes, we have adopted another 
approximation which might breakdown at asymptotic times. We have thus found that the 
fermion production is in general mildly suppressed by the back-reaction (see Figs. [TT] and [T2l) . 
At later times (r > 10^ T), the back-reaction of the produced top has been neglected and, 
in the same framework of approximations, we have found that the production of Z° is more 
strongly suppressed by the back-reaction (see Figs. [T^ and [T^ . 
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The study of the system where both fermion and boson backreaction effects are included 
has shown a slower boson production, which takes over almost an order of magnitude later 
in time, against the naive expectation of some kind of factorization in the production events. 
We however remark that we have chosen a value of the Higgs mass Mh = 500 GeV which 
corresponds to a large production of the top in order to study its effect on the production of 
the Z". For a more realistic (presumably smaller) value of M^, one should instead consider 
the quarks and leptons with the higher production rates (as shown in Table [T]). 

Future investigations, apart from being focused on more appropriate values of the Higgs 
mass as, for example, an awaited discovery at LHC would provide, should address aspects 
which we have found beyond our possibility because of the approximation schemes adopted. 
In particular, one should investigate the production and dissipative dynamics in the late 
(asymptotic) times of the evolution and evaluate the rescattering phenomena associated 
with the produced particles. Keeping in mind that we are mainly considering a possible 
"starting" regime of very small Higgs oscillations, the rescattering could lead to a thermal 
background characterized by a very low temperature. 

We conclude by noting that, although suppressed when the back-reaction is properly 
included in the analysis, particle production by parametric resonance with an oscillating 
Higgs field remains a remarkable effect with phenomenological relevance and could be an 
(indirect) way of testing the time-dependence of the constants of the Standard Model. 
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